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A theoretical approach to thermal noise caused by an inhomogeneously distributed 
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— Physical insight by the advanced modal expansion 
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We modified the modal expansion, which is the traditional method used to calculate thermal noise. 
This advanced modal expansion provides physical insight about the discrepancy between the actual 
thermal noise caused by inhomogeneously distributed loss and the traditional modal expansion. This 
discrepancy comes from correlations between the thermal fluctuations of the resonant modes. The 
thermal noise spectra estimated by the advanced modal expansion are consistent with the results of 
measurements of thermal fluctuations caused by inhomogeneous losses. 
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I. INTRODUCTION 

Thermal fluctuation is one of the fundamental noise 
sources in precise measurements. For example, the sen- 
sitivity of interferometric gravitational wave detectors 
d S B 0] is limited by the thermal noise of the mechan- 
ical components. The calculated thermal fluctuations of 
rigid cavities have coincided with the highest laser fre- 
quency stabilization results ever obtained [H, It is 
important to evaluate the thermal motion for studying 
the noise property. The (traditional) modal expansion 
has been commonly used to calculate the thermal 
noise of elastic systems. However, recent experiments 
d, H, EH EH E2, EH have revealed that modal expan- 
sion is not correct when the mechanical dissipation is 
distributed inhomogeneously. In some theoretical stud- 
ies [3, El, EH E3], calculation methods that are com- 
pletely different from modal expansion have been devel- 
oped. These methods are supported by the experimental 
results of inhomogeneous loss d, EH EH EH- However, 
even when these method were used, the physics of the 
discrepancy between the actual thermal noise and the 
traditional modal expansion was not fully understood. 

In this paper, another method to calculate the thermal 
noise is introduced [17| . This method, advanced modal 
expansion, is a modification of the traditional modal ex- 
pansion (this improvement is a general extension of a 
discussion in Ref. [18]). The thermal noise spectra es- 
timated by this method are consistent with the results 
of experiments concerning inhomogeneous loss [1, E3| ■ It 
provides information about the disagreement between the 
thermal noise and the traditional modal expansion. We 
present the details of these topics in the following sec- 
tions. 
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II. OUTLINE OF ADVANCED MODAL 
EXPANSION 

A. Review of the traditional modal expansion 

The thermal fluctuation of the observed coordinate, 
X, of a linear mechanical system is derived from the 
fluctuation-dissipation theorem [l9|, H3, HH , 



Gx(f) = -^lm[H x (oj) 



H x (lu) = 
X(co) = 



^ poo 

2^ 



X(t) cxp(—iujt)dt, 



(1) 
(2) 
(3) 



where f(— lu/2tt), t, fee and T, are the frequency, time, 
Boltzmann constant and temperature, respectively. The 
functions (Gx, Hx, and F) are the (single-sided) power 
spectrum density of the thermal fluctuation of X, the 
transfer function, and the generalized force, which corre- 
sponds to X. In the traditional modal expansion [7j], in 
order to evaluate this transfer function, the equation of 
motion of the mechanical system without any loss is de- 
composed into those of the resonant modes. The details 
are as follows: 

The definition of the observed coordinate, X, is de- 
scribed as 



X(t) = J u(r,t) • P(r)dS, 



(4) 



where u is the displacement of the system and P is a 
weighting function that describes where the displacement 
is measured. For example, when mirror motion is ob- 
served using a Michelson interferometer, as in Fig. [U X 
and u represent the interferometer output and the dis- 
placement of the mirror surface, respectively. The vector 
P is parallel to the beam axis. Its norm is the beam- 
intensity profile EH • The equation of motion of the me- 
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FIG. 1: Example of the definition of the observed coordinate, 
X, in Eq. @. The mirror motion is observed using a Michel- 
son interferometer. The coordinate X is the output of the 
interferometer. The vector u represents the displacement of 
the mirror surface. The field P is parallel to the beam axis. 
Its norm is the beam-intensity profile [TJ|]. 



chanical system without dissipation is expressed as 
fp'n 

p—-C[u]=F{t)P(r) 1 



(5) 



where p is the density and C is a linear operator. The 
first and second terms on the left-hand side of Eq. ([5]) 
represent the inertia and the restoring force of the small 
elements in the mechanical oscillator, respectively. The 
solution of Eq. §5§ is the superposition of the basis func- 
tions, 



u(r,t) = y^n> n (r)g n (f). 



(6) 



The functions, w n and q n , represent the displacement 
and time development of the n-th resonant mode, re- 
spectively. The basis functions, w n , are solutions of the 
eigenvalue problem, written as 



C[w n (r)} = -puj n 2 w n (r), 



(7) 



where oj n is the angular resonant frequency of the n-th 
mode. The displacement, w n , is the component of an 
orthogonal complete system, and is normalized to satisfy 
the condition 



w n (r) ■ P(r)dS = 1. 



(8) 



J pw n (r) ■ w k (r)dV 



(9) 



The parameter m n is called the effective mass of the 
mode H HH, [23|, [24I ■ The tensor S n k is the Kronecker's 
(5-symbol. 

Putting Eq. © into Eq. (J), we obtain a relationship 
between X and q n using Eq. (|8|), 



*(*) = 5>„(t). 



(10) 



In short, coordinate X is a superposition of those of the 
modes, q n . In order to decompose the equation of mo- 
tion, Eq. J5]), Eq. ([6|) is substituted for u in Eq. ((5]). 
Equation ([5]) is multiplied by w n and then integrated 
over all of the volume using Eqs. (©, © and ©. The 
result is that the equation of motion of the n-th mode, 
q n , is the same as that of a harmonic oscillator on which 
force F(t) is applied. After modal decomposition, the 
dissipation term is added to the equation of each mode. 
The equation of the n-th mode is written as 



m n uj q n + m n uj n [1 + i(j) n {uj)]q n 



F, 



(11) 



in the frequency domain. The function <f) n is the loss 
angle, which represents the dissipation of the n-th mode 
[21 . The transfer function, Hx, derived from Eqs. ([2]), 
(fTOf and (fTTj) is the summation of those of the modes, 



E 



-m„w 2 + m n u) n 2 [l + i</>„(w)] 



■ (12) 



According to Eqs. fl} and (fTSj) . the power spectrum 
density, Gx, is the summation of the power spectrum, 
Gq n , of q n , 



n 

Ak B T 



E 



m n uj (uj 2 - uj n 2 ) 2 + uj n 4 (j) n 2 (u>) 



(13) 



B. 



Equation of motion in an advanced modal 
expansion 



In the traditional modal expansion, the dissipation 
term is introduced after decomposition of the equation 
of motion without any loss. On the contrary, in an ad- 
vanced modal expansion, the equation with the loss is 
decomposed [T3, HE] . If the loss is sufficiently small, the 
expansion process is similar to that in the perturbation 
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theory of quantum mechanics [26j . The equation of q n is 
expressed as 

- m n uj 2 q n + m„cj„ 2 [l + i(f> n (u>)]q n 



^ ia„ fc (w)q fc = F, (14) 
fc^n 



(15) 



The third term in Eq. (fl4"l) is the difference between 
the advanced, Eq. (fT4"|) , and traditional, Eq. (|11[) . modal 
expansions. Since this term is a linear combination of the 
motions of the other modes, it represents the couplings 
between the modes. The magnitude of the coupling, a. n u, 
depends on the property and the distribution of the loss 
(described below). 



C. Details of coupling 

Let us consider the formulae of the couplings caused by 
the typical inhomogeneous losses, the origins of which ex- 
ist outside and inside the material (viscous damping and 
structure damping, respectively) [17| . Regarding most of 
the external losses, for example, the eddy-current damp- 
ing and residual gas damping are of the viscous type [7| . 
The friction force of this damping is proportional to the 
velocity. Inhomogeneous viscous damping introduces a 
friction force, iujpT(r)u(r), into the left-hand side of the 
equation of motion, Eq. ([5]), in the frequency domain. 
The function T(> 0) represents the strength of the damp- 
ing. The equation of motion with the dissipation term, 
iuipT(r)u(r), is decomposed. Since the loss is small, the 
basis functions of the equation without loss are available 
[26|. Equation ([6]) is put into the equation of motion 
along with the inhomogeneous viscous damping. This 
equation multiplied by w n is integrated. The coupling of 
this dissipation is written in the form 



a„ fc = oj / pT(r)w n (r) ■ w k {r)dV = a kT: 



(16) 



In most cases, the internal loss in the material is ex- 
pressed using the phase lag, </>(> 0), between the strain 
and the stress 0- The magnitude of the dissipation is 
proportional to this lag. The phase lag is almost con- 
stant against the frequency 0] in many kinds of mate- 
rials (structure damping). In the frequency domain, the 
relationship between the strain and the stress (the gener- 
alized Hooke's law) in an isotropic elastic body is written 

as Bin nana 



E [l + ^(r)} 



1 



1 



2 °r 



uuSi, 



(rK-, 



(17) 
(18) 



= [1- 

1 / dui duj 

13 2 \dxj dxi 
where E Q is Young's modulus and a is the Poisson ratio; 
Gij and Uij are the stress and strain tensors, respectively. 
The tensor, 5^, is the real part of the stress, It rep- 
resents the stress when the structure damping vanishes. 
The value, m, is the z-th component of u. The equation 
of motion of an elastic body [27| in the frequency domain 
is expressed as 



pu 2 Ui 



(19) 



where Pi is the i-th component of P. From Eqs. (jT7|) 
and (|19p . an inhomogeneous structure damping term is 
obtained, —\'Y^jd(j){r)d' i j/dxj. The equation of motion 
with the inhomogeneous structure damping is decom- 
posed in the same manner as that of the inhomogeneous 
viscous damping. The coupling is calculated using inte- 
gration by parts and Gauss' theorem [27l |. 
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E 



dV = a kr 



(20) 



r 



where w n ,i and rij are the i-th components of «;„ and 
the normal unit vector on the surface. The tensors, w n ,ij 
and a n ,iji are the strain and stress tensors of the n-th 
mode, respectively. In order to calculate these tensors, 
w n .i is substituted for m in Eqs. (jTTJ) and (fT5|) with 
= 0. Equation (j2TJ|) is valid when the integral of the 
function, ■ w n ^4><Jk,ijnj , on the surface of the elastic 
body vanishes. For example, the surface is fixed (u>„,i = 



0) or free (V 0*,- 



0) (27|. 



The equation of motion in the advanced modal expan- 
sion coincides with that in the traditional modal expan- 
sion when all of the couplings vanish. A comparison be- 
tween Eqs. and (HHJ) shows that in viscous damping 
all a nk (n ^ k) are zero when the dissipation strength, 
T(r), does not depend on the position, r. In the case of 
structure damping, from Eqs. ©, ©, Q and (fP5)) . the 



stress, a[j, without dissipation satisfies 



E 

3 



dxj 



(21) 



According to Eq. J9J), Eq. (f2T[) is decomposed without 
any couplings. From Eq. (|21[) and the structure damping 
term, — i XZj ^4'( r )^ij I dxj, the conclusion is derived; all 
of the couplings in the structure damping vanish when 
the loss amplitude, 4>, is homogeneous. In summary, the 
inhomogeneous viscous and structure dampings produce 
mode couplings and destroy the traditional modal expan- 



does not change while the resonant motion decays. On 
the other hand, if the dissipation is inhomogeneous, the 
motion near the concentrated loss decays more rapidly 
than the other parts. The shape of the displacement be- 
comes different from that of the original resonant mode. 
This implies that the other modes are excited, i.e. the 
energy of the original mode is leaked to the other modes. 
This energy leakage represents the couplings in the equa- 
tion of motion. 

It must be noticed that some kinds of " homogeneous" 
loss cause the couplings. For example, in thermoelastic 
damping [H, [29|, [30, HU , which is a kind of internal loss, 



the energy components of the shear strains, 



are not dissipated. The couplings, a nk , do not have any 
terms that consist of the shear strain tensors. The cou- 
pling formula of the homogeneous thermoelastic damp- 
ing is different from Eq. (j2"0)) with the constant <f>. The 
couplings are not generally zero, even if the thermoelas- 
tic damping is uniform. The advanced, not traditional, 
modal expansion provides a correct evaluation of the "ho- 
mogeneous" thermoelastic damping. In this paper, how- 
ever, only coupling caused by inhomogeneous loss is dis- 
cussed. 



D. Thermal- noise formula of advanced modal 
expansion 



The reason why the inhomogeneity of the loss causes 
the couplings is as follows. Let us consider the decay mo- 
tion after only one resonant mode is excited. If the loss 
is uniform, the shape of the displacement of the system 



In the advanced modal expansion, the transfer func- 
tion, H x , is derived from Eqs. {2J, pH]) . and fli]) (since 
the dissipation is small, only the first-order of a nk is con- 
sidered [32l|). 



J 



i 



-m n u) 2 + m n uj n 2 (l + i(j) n ) 



E 

k^n 



Milt 



-m n uj 2 + m n LVn 2 (l + i4> n )][-m k Lu 2 + m k ui k 2 {l + i(j) k )] 



(22) 
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Putting Eq. into Eq. flU, the formula for the ther- 
mal noise is obtained. In the off-resonance region, where 
| — uj 2 + uj n 2 \ 3> w„ 2 n (o;) for all n, this formula approx- 
imates the expression 



Gx(f) = 



4fc B r u; n 2 (i> n (uj) 

^ m n uj (uj 2 - U! n 2 ) 2 



1 1 



4fc B T 



f-f m n m k uj (uj 2 



Un 2 )(u 2 - UJk 2 ) 



(23) 



The first term is the same as the formula of the tradi- 
tional modal expansion, Eq. (|13p . 

The interpretation of Eq. (|23[) is as follows. The power 
spectrum density of the thermal fluctuation force of the 
n-th mode, Gf„ , and the cross-spectrum density between 
F n and Fu, GF n F k , are evaluated from Eq. (fT4")) and the 
fluctuation-dissipation theorem 0, [2l[ , 



G Fn (f) = 4fc B T 
G FnFk (f) = 4fc B T 



m n uj n 2 (/) n (uj) 

UJ 

UJ 



(24) 
(25) 



The power spectrum density, Gf„ , is independent of a n k- 
On the other hand, Gp n F k depends on a n k- Having the 
correlations between the fluctuation forces of the modes, 
correlations between the motion of the modes must also 
exist. The power spectrum density of the fluctuation 
of q n , G qn , and the cross-spectrum density between the 



fluctuations of q n and qt, G qnqk , are described as 20, 21 



The calculation used to evaluate the mean square is writ- 
ten as 



21] 



X 2 



Gx(f)df 



1 f°° 

— I G x (u)duj 

47r J-oo 

k B T f°° lm[H x (uj)} 



du> 



7T J-oo & 

k B TRe[H x (0)}. 



(29) 



Since the transfer function, H x , is the ratio of the Fourier 
components of the real functions, the value Hx(0) is a 
real number. The functions <j> n and a n k, which cause 
the imaginary part of Hx, must vanish when uj is zero 
[2l| . The correlations do not affect the mean square of 
the thermal fluctuation. Equation (|29p is rewritten using 
Eq. (H2J) as 



X 2 



E 



k B T 



m n uj n 



2 ' 



(30) 



Equation (I30|) is equivalent to the prediction of the 
equipartition principle. 

The calculation of the formula of the advanced modal 
expansion, Eq. (1231) , is more troublesome than that of the 
other methods I14L Il5l . , , which are completely dif- 
ferent from the modal expansion, when many modes con- 
tribute to the thermal motion. However, the advanced 
modal expansion gives clear physical insight about the 
discrepancy between the thermal motion and the tradi- 
tional modal expansion, as shown in Sec. HVl It is diffi- 
cult to find this insight using other methods. 



G qn (f) 

Gq n q k (/) 



4fc B r UJ n 2 (f) n (uj) 

m n uj (uj 2 - uj n 2 ) 2 ' 

4fc B T Onk 

m n mkUJ (uj 2 - uj n 2 )(uj 2 - UJk 2 ) 



(26) 
(27) 



under the same approximation of Eq. (|23p . The first and 
second terms in Eq. (|23|) are summations of the fluctu- 
ation motion of each mode, Eq. (j2"6")l , and the correla- 
tions, Eq. (|27p. respectively. In conclusion, inhomoge- 
neous mechanical dissipation causes mode couplings and 
correlations of the thermal motion between the modes. 

In order to check wheather the formula of the thermal 
motion in the advanced modal expansion is consistent 
with the the equipartition principle, the mean square of 
the thermal fluctuation, X 2 , which is an integral of the 
power spectrum density over the whole frequency region, 
is evaluated. This mean square is derived from Eq. ([T]) 
using the Kramers-Kronig relation [U, |33| , 



Re[H x (uj)} 



MHx(0\ 



UJ 



(28) 



III. EXPERIMENTAL CHECK 

In order to test the advanced modal expansion exper- 
imentally, our previous experimental results concerning 
oscillators with inhomogeneous losses 0, [UJ are com- 
pared with an evaluation of the advanced modal expan- 
sion [l7l |. In an experiment involving a drum (a hol- 
low cylinder made from aluminum alloy as the prototype 
of the mirror in the interferometer) with inhomogeneous 
eddy-current damping by magnets [12j , the measured val- 
ues agreed with the formula of the direct approach [3] , 
Eq. (6) in Ref. d2j. This expression is the same as that 
of the advanced modal expansion [171 ]. 

FigureOpresents the measured spectra of an aluminum 
alloy leaf spring with inhomogeneous eddy-current damp- 
ing [8j . The position of the magnets for the eddy-current 
damping and the observation point are indicated above 
each graph. In the figures above each graph, the left 
side of the leaf spring is fixed. The right side is free. 
The open circles in the graphs represent the power spec- 
tra of the thermal motion derived from the measured 
transfer functions using the fluctuation-dissipation theo- 
rem. These values coincide with the directly measured 
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FIG. 2: Comparison between the estimated advanced modal expansion and the experimental results of an aluminum alloy 
leaf spring with inhomogeneous eddy-current damping Q. The position of the magnets for the eddy-current damping and the 
observation point are indicated above each graph. In the figures above each graph, the left side of the leaf spring is fixed. The 
right side is free. The open circles in the graphs represent the power spectra of the thermal motion derived from the measured 
transfer functions using the fluctuation-dissipation theorem. These values coincide with the directly measured thermal-motion 
spectra ,8j. The solid lines are estimations using the advanced modal expansion. As a reference, an evaluation of the traditional 
modal expansion is also given (dashed lines). 



7 



thermal- motion spectra [8|. The solid lines are estima- 
tions using the advanced modal expansion (the correla- 
tions derived from Eqs. (fTBf and ([27]) are almost perfect 
[l7l]). As a reference, an evaluation of the traditional 
modal expansion is also given (dashed lines). The re- 
sults of a leaf-spring experiment are consistent with the 
advanced modal expansion. Therefore, our two experi- 
ments support the advanced modal expansion. 



IV. PHYSICAL INSIGHT GIVEN BY THE 
ADVANCED MODAL EXPANSION 

The advanced modal expansion provides physical in- 
sight about the disagreement between the real thermal 
motion and the traditional modal expansion. Here, let 
us discuss the three factors that affect this discrepancy: 
the number of the modes, the absolute value and the sign 
of the correlation. 



A. Number of modes 

Since the difference between the advanced and tradi- 
tional modal expansions is the correlations between the 
multiple modes, the number of the modes affects the mag- 
nitude of the discrepancy. If the thermal fluctuation is 
dominated by the contribution of only one mode, this 
difference is negligible, even when there are strong corre- 
lations. On the other hand, if the thermal motion con- 
sists of many modes, the difference is larger when the 
correlations are stronger. 

Examples of the one-mode oscillator are given in Fig. 
[2J The measured thermal motion spectra of the leaf 
spring with inhomogeneous losses below 100 Hz were the 
same as the estimated values of the "traditional" modal 
expansion. This is because these fluctuations were dom- 
inated by only the first mode (about 60 Hz) . As another 
example, let us consider a single-stage suspension for a 
mirror in an interferometric gravitational-wave detector. 
The sensitivity of the interferometer is limited by the 
thermal noise of the suspensions between 10 Hz and 100 
Hz. Since, in this frequency region, this thermal noise is 
dominated by only the pendulum mode 0], the thermal 
noise generated by the inhomogeneous loss agrees with 
the traditional modal expansion. It must be noticed that 
the above discussion is valid only when the other sus- 
pension modes are negligible. For example, when the 
laser beam spot on the mirror surface is shifted, the two 
modes (pendulum mode and mirror rotation mode) must 
be taken into account. In such cases, the inhomogeneous 
loss causes a disagreement between the real thermal noise 
of the single-stage suspension and the traditional modal 
expansion [34j . 

The discrepancy between the actual thermal motion 
and the traditional modal expansion in the elastic modes 
of the mirror [24[ is larger than that of the drum, the 
prototype of the real mirror in our previous experiment 



[121 ]. One of the reasons is that the thermal motion of 
the mirror (rigid cylinder) consists of many modes [22I 

tThe drum (hollow cylinder) had only two modes 
. Since the number of modes that contribute to the 
thermal noise of the mirror in the interferometer increases 
when the laser beam radius becomes smaller [22|, [Hj], 
the discrepancy is larger with a narrower beam. This 
consideration is consistent with our previous calculation 

B. Absolute value of the correlation 

In Eq. ([27)) , the absolute value of the cross-spectrum 
density, G 9 „ (?fc , is proportional to that of the coupling, 
a n k ■ Equations (fll>)) and ([2H)l show that the coupling de- 
pends on the scale of the dissipation distribution. A sim- 
ple example of viscous damping is shown in Fig. [3] Let 
us consider the absolute value of a n k when the viscous 
damping is concentrated (at around r v ; s ) in a smaller 
volume (AV) than the wavelengths of the n-th and fc-th 
modes. An example of this case is (A) in Fig. [31 It 
is assumed that the vector w n (r V i S ) is nearly parallel to 
Wk(r v i s ). The absolute value of the coupling is derived 
from Eqs. (fT5)) and (fTB) as 

vis jl^n \T vis J ' ^k\^vis )AV| 

~ (r v is ) I w n (r vis )\ 2 AV 

x ^pr(r v i s )|™ fc (r vis )| 2 Ay 

~ \fa nn a kk = vWAniwS' (31) 

The absolute value of the cross-spectrum is derived from 
Eqs. (EH), J37J, and as 

\Gq n q k \ ~ \/ G qn G qk . (32) 

In short, the correlation is almost perfect [35[. On the 
other hand, if the loss is distributed more broadly than 
the wavelengths, the coupling, i.e. the correlation, is 
about zero, 

|G g „ g J~0. (33) 

The dissipation in the case where the size is larger than 
the wavelengths is equivalent to the homogeneous loss. 
An example of this case is (B) in Fig. [3] Although 
the above discussion is for the case of viscous damping, 
the conclusion is also valid for other kinds of dissipa- 
tion. When the loss is localized in a small region, the 
correlations among many modes are strong. The loss in 
a narrower volume causes a larger discrepancy between 
the actual thermal motion and the traditional modal ex- 
pansion. This conclusion coincides with our previous cal- 
culation of a mirror with inhomogeneous loss [24( . 

C. Sign of correlation 

The sign of the correlation depends on the frequency, 
the loss distribution, and the position of the observation 
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FIG. 3: Example for considering the absolute value of the 
coupling. There are the n-th and fc-th modes, w n and Wk, of 
a bar with both free ends. The vertical axis is the displace- 
ment. The dashed horizontal lines show the bar that does not 
vibrate. When only the grey part (A), which is narrower than 
the wavelengths on the left-hand side, has viscous damping, 
the absolute value of the coupling, Eq. (|16|l . is large. Because 
the signs of w n and Wk do not change in this region. If vis- 
cous damping exits only in the hatching part (B), which is 
wider than the wavelengths on the right-hand side, the cou- 
pling is about zero, because, in this wide region, the sign of 
the integrated function in Eq. (|16p . which is proportional to 
the product of w n and Wk, changes. 



area. The position dependence provides a solution to the 
inverse problem: an evaluation of the distribution and 
frequency dependence of the loss from measurements of 
the thermal motion. 



1. Frequency dependence 

According to Eq. the sign of the correlation re- 

verses at the resonant frequencies. For example, in calcu- 



lating the double pendulum [18J, experiments involving 
the drum [l^] and a resonant gravitational wave detector 
with optomechanical readout [l(| , this change of the sign 
was found. In some cases, the thermal-fluctuation spec- 
trum changes drastically around the resonant frequen- 
cies. A careful evaluation is necessary when the obser- 
vation band includes the resonant frequencies. Examples 
are when using wide-band resonant gravitational- wave 
detector s f3rj [33, [H, Hi| , and thermal-noise interferome- 
ters [ll|, Il3|. The reason for the reverse at the resonance 
is that the sign of the transfer function of the mode with 
a small loss from the force (F n ) to the motion (q n ), H n 
in Eq. (fl2|) [cx (— uj 2 + w„ 2 ) _1 ], below the resonance is 
opposite to that above it. 

Since the sign of the correlation changes at the reso- 
nant frequencies, the cross-spectrum densities, the second 
term of Eq. (|23|) . make no contribution to the integral 
of the power spectrum density over the whole frequency 
region, i.e. the mean square of the thermal fluctuation, 
X 2 , as shown in Sec. Ill Dl Therefore, the consideration 
in Sec. Ill Dl indicates that a reverse of the sign of the 
correlation conserves the equipartition principle, a fun- 
damental principle in statistical mechanics. 



2. Loss and observation area position dependence 

According to Eqs. (flB]) and (j2"0)) . and the normaliza- 
tion condition, Eq. (JS]) [4(j, the sign of the coupling, 
a n k, depends on the loss distribution and the position 
of the observation area. A simple example is shown in 
Fig. |U Owing to this normalization condition, near the 
observation area, the basis functions, w n , are similar in 
most cases. On the contrary, in a volume far from the 
observation area, w n is different from each other in many 
cases. From Eqs. ©, (HHJ), & and ([25]), when the loss 
is concentrated near to the observation area, most of the 
couplings (and the correlations between the fluctuation 
forces of the modes, Gp n F h ) are positive. On the other 
hand, when the loss is localized far from the observation 
area, the numbers of the negative couplings and Gp n F k 
are about the same as the positive ones. In such a case, 
most of the couplings between the n-th and (n ± l)-th 
modes (and Gf„p n±1 ) are negative. These are because 
the localized loss tends to apply to the fluctuation force 
on all of the modes to the same direction around itself. 
Equation (|2~T)) indicates that the sign of the correlation, 
Gq nqk , is the same as that of the coupling, a n k, below 
the first resonance. In this frequency band, the thermal 
motion is larger and smaller than the evaluation of the 
traditional modal expansion if the dissipation is near and 
far from the observation area, respectively. This conclu- 
sion is consistent with the qualitative discussion of Levin 
[Tl |. our previous calculation of the mirror and the 
drum experiment fl2| . 
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Resonant modes of a bar 
with both free ends 

Observation 
point 

1) Bar 




1st mode 



2nd mode 



3rd mode 



FIG. 4: Example for considering the sign of the coupling. 
There are the lowest three modes, w n , of a bar with both 
free ends. The vertical axis is the displacement. The dashed 
horizontal lines show the bar that does not vibrate. The ob- 
servation point is at the right-hand side end. The normal- 
ization condition is Eq. (|8]) [40| . The sign and shape of the 
displacement of all the modes around the observation point 
are positive and similar, respectively. On the contrary, at the 
left-hand side end, the sign and shape of the n-th mode are 
different from each other in many cases. From Eqs. (JSJ), (JTSJ, 
(|20l) . when the loss is concentrated near to the observation 
area, most of the couplings are positive. On the other hand, 
when the loss is localized far from the observation area, the 
number of the negative couplings is about the same as the 
positive one. In such a case, most of the couplings between 
the n-th and (n ± l)-th modes are negative. 



3. Inverse problem 

The above consideration about the sign of the coupling 
gives a clue to solving the inverse problem: estimations 
of the distribution and frequency dependences of the loss 
from the measurement of the thermal motion. Since the 
sign of the coupling depends on the position of the ob- 
servation area and the loss distribution, a measurement 
of the thermal vibrations at multiple points provides in- 
formation about the couplings, i.e. the loss distribution. 
Moreover, multiple-point measurements reveal the loss 
frequency dependence. Even if the loss is uniform, the 
difference between the actual thermal motion and the 
traditional modal expansion exits when the expected fre- 
quency dependence of the loss angles, 0„(w), is not cor- 
rect [18J. The measurement at the multiple points shows 
whether the observed difference is due to an inhomoge- 
neous loss or an invalid loss angle. This is because the 
sign of the difference is independent of the position of the 
observation area if the expected loss angles are not valid. 

As an example, our leaf-spring experiment Q is dis- 
cussed. The two graphs on the right (or left) side of Fig. 
[5] show thermal fluctuations at different positions in the 
same mechanical system. The spectrum is smaller than 
the traditional modal expansion. The other one is larger. 
Thus, the disagreement in the leaf-spring experiment was 
due to inhomogeneous loss, not invalid loss angles. When 
the power spectrum had a dip between the first (60 Hz) 
and second (360 Hz) modes, the sign of the correlation, 
G qiQ2 , was negative. According to Eq. P7|) . the sign 
of the coupling, a\2 1 was positive. The loss was concen- 
trated near to the observation point when a spectrum dip 
was found. The above conclusion agrees with the actual 
loss shown in Fig. O 



V. CONCLUSION 

The traditional modal expansion has frequently been 
used to evaluate the thermal noise of mechanical systems 
0. However, recent experimental research [H, |T(| El, 
H2I . [IH has proved that this method is invalid when the 
mechanical dissipation is distributed inhomogeneously. 
In this paper, we introduced a modification of the modal 
expansion [l7l . Il8j . According to this method (the ad- 
vanced modal expansion), inhomogeneous loss causes cor- 
relations between the thermal fluctuations of the modes. 
The fault of the traditional modal expansion is that these 
correlations are not taken into account. Our previous 
experiments concerning the thermal noise of the 

inhomogeneous loss support the advanced modal expan- 
sion. 

The advanced modal expansion gives interesting phys- 
ical insight about the difference between the actual ther- 
mal noise and the traditional modal expansion. When 
the thermal noise consists of the contributions of many 
modes, the loss is localized in a narrower area, which 
makes a larger difference. When the thermal noise is 
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dominated by only one mode, this difference is small, 
even if the loss is extremely inhomogeneous. The sign of 
this difference depends on the frequency, the distribution 
of the loss, and the position of the observation area. It 
is possible to derive the distribution and frequency de- 
pendence of the loss from measurements of the thermal 
vibrations at multiple points. 

There were many problems concerning the thermal 
noise caused by inhomogeneous loss. Our previous work 
[li _L_i and this research solved almost all of these 
problems: a modification of the traditional estimation 
method (in this paper), experimental checks of the new 
and traditional estimation methods and a confirmation of 
the new methods ([__, ___} and this paper), an evaluation 



of the thermal noise of the gravitational wave detector 
using the new method and a consideration of the 
physical properties of the discrepancy between the actual 
thermal noise and the traditional estimation method (in 
this paper). 
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